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. Abstract 

o ■ 

■~«.^ ■ We consider the mechanical equilibrium and stability of matter trapped 

i-C ' in the magnetosphere of a rapidly rotating star. Assuming a dipolar mag- 

netic field and arbitrary inclination of the magnetic axis with respect to 
^ ■ the axis of rotation we find stable equilibrium positions a) in a (warped) 

5-H ' disk roughly aligned with the magnetic equatorial plane and b) at two lo- 

, cations above and below the disk, whose distance from the star increases 

' with decreasing inclination angle between dipole and rotation axis. The 

distribution of matter is not strongly affected by allowing for a spatial off- 
set of the magnetic dipole. These results provide a possible explanation 
' for some observations of corotating localized mass concentrations in hot 

' magnetic stars. 

1 Introduction 



Spectroscopic observations indicate the presence of corotating circumstellar mat- 
ter around chemically peculiar stars (hehum-strong, helium- weak, and Ap/Bp 
stars) with strong magnetic fields (e.g.. Shore & Adelman 1981, Walborn 1982, 
Shore et al. 1990, Groote & Hunger 1997, Smith et al. 1998, Smith 2001). 
Such accumulation of mass often coexists with magnetically confined outflows 
(Barker et al. 1982, Brown et al. 1985, Shore et al. 1987, Smith & Groote 
2001). Shore and Brown (1990, see also Shore 1987) suggested that the corotat- 
ing plasma is trapped in a magnetically closed region near the magnetic equator 
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while the stellar wind forms jet-like outflows along open fleld lines around the 
magnetic poles. This is similar to the earlier concept of Mestel (1968, see also 
Mestel & Spruit 1987) of 'dead zones' and 'wind zones' and also reminiscent 
to the structure of the solar corona during activity minimum (e.g., Pneuman 
& Kopp 1971). Recently, Ud-Doula & Owocki (2002) have carried out two- 
dimensional MHD simulations illustrating the development of open and closed 
field topologies around a non-rotating hot star with a line-driven wind. 

A mechanism for the accumulation of circumstellar matter in equatorial 
disks around Be stars has been proposed by Bjorkman & Cassinelli (1993, see 
also Cassinelli et al. 2002, and references therein). They suggest that the strong 
magnetic field forces the radiatively driven outflowing matter to follow the dipo- 
lar fleld lines and collide near the magnetic equator. Since the colliding flows 
are supersonic, shock heating leads to UV and X-ray emission. After the mate- 
rial has cooled sufficiently, it accumulates near the magnetic equator and forms 
a disk. Babel and Montmerle (1997) have proposed a similar mechanism for 
corotating circumstellar matter and X-ray emission from Ap/Bp stars. For the 
matter to become magnetically trapped in a closed part of the magnetosphere, 
the kinetic energy density of the wind has to be smaller than the magnetic en- 
ergy density. Depending on the strength of the magnetic field, wind speed and 
mass loss rate, estimates yield an extension of the corotating magnetosphere of 
up to 10 stellar radii. 

Most of the previous theoretical considerations of the equilibrium and dy- 
namics of circumstellar matter around magnetic stars considered aligned rota- 
tors, for which the magnetic axis and the rotation axis coincide. In this case, 
the trapped matter accumulates near the equatorial plane. The distribution 
of the magnetospheric matter in an oblique rotator with arbitrary relative in- 
clination of the two axes is far less clear. Here we consider this general case 
by studying the force equilibrium and the stability of circumstellar matter in 
the closed magnetosphere of an oblique rotator. Assuming that the magnetic 
energy density is much larger than the kinetic energy density of rotation, which, 
in turn, is much larger than the thermal energy density of the plasma, we can 
assume the magnetic field to be fixed and ignore the thermal pressure. The 
equilibrium problem is then reduced to finding the locations where the (vector) 
sum of the gravitational and the centrifugal force is perpendicular to a given 
magnetic field line. A similar approach has been taken by Nakajima (1985), who 
considered the minima of the effective potential energy along a given magnetic 
fleld line. Nakajima apparently determined the absolute minima of the poten- 
tial and thus could only find one equilibrium locus per field line. In our work, 
however, we show that there are more locations where the trapped plasma can 
accumulate. We thus are able to predict the complete spatial distribution of the 
magnetospheric plasma in a rapidly rotating magnetic star. 

The paper is organized as follows. Sec. |21 describes the method to determine 
the force equilibrium and its stability. Sec. |31 gives results for the special cases 
of aligned and perpendicular rotators as well as for other relative inclinations 
of the rotation and magnetic axes. We also consider offset dipoles. We discuss 
the results in Sec. 0] and present our conclusions in Sec. El 
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2 Force equilibrium and stability 



We consider the closed part of the magnetosphere of a rapidly rotating star, i.e., 
the region inside the Alfven radius defined by the equality of stellar wind speed 
and Alfven velocity. We assume the scaling 



« vl « vl , (1) 



where fg, uo, and v\ are the local sound speed, rotational speed (in an inertial 
frame), and Alfven speed, respectively. This scaling corresponds to a situation 
where the magnetic energy density is much larger than the kinetic energy den- 
sity of the rotational motion, which, in turn, is much larger than the thermal 
energy density of the magnetospheric matter. Under these conditions, which are 
relevant for the observed stellar magnetospheres out to about 10 stellar radii 
(Nakajima 1985; Shore 1987; Babel & Montmerle 1997), the magnetic field is 
largely unaffected by the presence of circumstellar matter. In the limit of very 
large electrical conductivity, the gas is attached to the magnetic field lines and 
the distribution of the circumstellar matter along a given field line is determined 
by the gravitational and the centrifugal forces (in a corotating frame of refer- 
ence). Since under the scaling given by Eq. (QJ the field lines are practically 
rigid (a minimal distortion of the field can balance any force perpendicular to 
the field lines), the force equilibrium for a mass element 5m is given by the 
balance of the components of the gravitational and centrifugal forces tangential 
to the local field line, 

(Fc + Fg) ■ B = , (2) 
where B is the magnetic field vector, 

Fc = 5m [e^ - eo (e^ • en)] (3) 

is the centrifugal force and 

M5m 

= - G—^ er (4) 

is the gravitational force. Here e,. denotes the radial unit vector, pointing from 
the center of the star to the mass element located at a radial distance r. eji is 
the unit vector along the axis of rotation, G denotes the gravitational constant 
and M the mass of the star. We consider a magnetic dipole field with a dipole 
moment, m, parallel to the unit vector em: 

B = ^^"g [3(er • m)er — m] . (5) 



Inserting Eqs. lO - © in Eq. yields 



1 1 f ^ 



- , (6) 
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which can be written as 



A {er ■ Gin) + (e^ ■ eo)(eo • e,„) (7) 

with 

^ = 2(l-^)-3(e„.eo)^ (8) 

The quantity GM/Vt^r^ represents the ratio of the gravitational to the cen- 
trifugal force. Kepler rotation would give a value of unity for this quantity. 

1/3 

(GM/SI^) = Tea is thus called the corotation radius. 

Eqs. 10) - I® show that the distribution of the equilibrium positions is scale- 
invariant for a fixed value of the angle ip between the magnetic and the rotational 
axes (eji -em = cosi/'). Therefore, writing r in units of rco allows us to obtain the 
result for any values of the stellar parameters Q. and M after proper rescaling 
with fco. 

To analyse the stability of a mass element located at an equilibrium position 
ro we consider the projection of the total force (F = Fc + Fq) on the tangent 
vector 1 of the field line for a small tangential displacement, 77, of a mass ele- 
ment from its equilibrium position along the corresponding magnetic field line. 
Denoting the equilibrium quantities by a subscript zero and neglecting second- 
and higher-order contributions in 77, we have 



F-l = (Fo + <5F) • (lo + (51) (9) 
= [Fo+r7(lo-V)F(ro)]-(lo + no77K-i), 

with 

F(ro + 77 lo) - Fo r? (Iq • V)F(ro) (10) 

and 

l(ro + 77lo) = lo + r;(lo- V)l(ro) (11) 
= lo + 77Ko^no. 

Ko is the local curvature radius and no is the normal vector to the field line 
through Tq. From Eq. @ we obtain 

F-l - {[(lo- V)F] -lo + Fo-noK^i}?] (12) 

to first order in 77, where we have used the equilibrium condition Fq Tq = 0. The 
neccessary and sufficient condition for a stable equilibrium is that the projection 
of F on 1 is directed towards the equilibrium position, viz. 

[(lo-V)F]-lo + Fo-noAto-^ <0. (13) 
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Figure 1: Hatched surfaces show the equihbrium positions for the case ip = 
(aUgned rotator) Both the magnetic dipole axis and the rotation axis are in 
z-direction. The left panel shows all equilibrium positions while the right panel 
shows only the stable ones. All lengths are given in units of the corotation 
radius. 

3 Results 

For the following illustrative examples we use parameter values corresponding 
to the AOp star IQ Aur (HD 34452, see Babel and Montmerle 1997), namely 
M = 4.8 Mq, = 5Rq and P,ot = 2.5 d ^ = 2.9 • lO^'^s"! and = 
13.4 i?0 — 2.6 i?*. Note that the results for cases with other values of M and 
Q can simply be found by rescaling with the corresponding value of , so that 
the results obtained here are, in general, valid for other stellar parameters as 
well. We first consider the special cases of aligned and perpendicular rotators 
and then show results for the oblique case. 

3.1 Aligned rotator: = 

From Eq. {Tj) we find for the case of parallel rotation and magnetic axes (em-eo = 
1) the equilibrium condition 




with the rotation axis directed along the z-axis. This equation has two solutions: 

1) cos 9 = Gr • efj = 0: this solution corresponds to the (coinciding magnetic 
and rotational) equatorial plane of the star. 




Figure 2: Same as Fig.n]for the perpendicular rotator (tp = 7r/2). The magnetic 
dipole is oriented along the x-axis. 



2) (e, • en) ^ 0: 

2/11 2 /rco\3 
COB ^ = 1-3 (-j 

The equilibrium locations described by this solution are chimney-shaped sur- 
faces above and below the equatorial plane. This solution only exists for r > 
(2/3)i/3rco. 

The left panel of Fig.nshows all equilibrium locations for the aligned rotator 
while the right panel gives only the stable equilibria. The 'chimney' equilibria 
are unstable while the equilibria in the equatorial plane are all stable outside 
the corotation radius. 

3.2 Perpendicular rotator: ip = 7i/2 

If the rotation axis and the magnetic axis are perpendicular to each other, the 
equilibrium condition, Eq. (UJ, becomes 



1 - 



3(e, • eo)' \ (e, • e„,) = . (15) 



Similar to the previous case, there are two solutions (see Fig.|21l: 

1) e,. • e,„ = 0: this solution corresponds to the equatorial plane with respect 
to the magnetic axis. 
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This corresponds to a chimney structure that is axisymmetric with respect to 
the rotation axis of the star. The solution only exists for r > Tco- 

For r ~ the co-latitude 9 is about 90°, corresponding to the base of 
the chimney at the equatorial plane. For increasing r, the co-latitude decreases 
down to a minimum value of 0niin — 35.3° (for r — > oo) so that at large distances 
from the star the chimneys form approximately coni with a constant opening 
angle 6'i„in- 

The stability analysis yields that both equilibrium surfaces are unstable in 
the close vicinity (out to 1-2 Tco) but are stable further out. 



2) • em 0: 

cos2 e = - 
3 



3.3 Oblique rotators 

If the inclination angle between the magnetic dipole and the rotation axis is 
oblique, we have to make use of the complete equilibrium condition, Eq. Q. 
Again we can find at least two independent solutions which, for small inclina- 
tions, describe in good approximation a disk-like and a chimney structure. For 
larger inclination angle, the shapes of the equilibrium surfaces are modified. 
This can be seen in Fig. 0] for = 0.257r, O.STStt, and 0.457r, respectively. For 
intermediate values, < ip < 0.5 tt, the disk becomes somewhat warped. In 
addition to the stable equilibrium points on the warped disk, there are also 
stable regions on the chimneys. Such regions exist for all values of ip, but they 
move away from the disk for decreasing inclination and reach infinity for = 0. 
The size and the curvature of these regions increases with increasing inclination. 
They are always completely separated from the disk, exept for the case -0 = 7r/2 
(Fig. 121) where the disk is directly connected to the stable regions on the chim- 
ney. In addition we find that the tilt angle 6 between the stable chimney regions 
and the rotation axis increases continously from ~ 0° for the aligned rotator up 
to 35.3° for the case of the perpendicular rotator. This implies that the new 
stable equilibrium regions can only exist within a specific angular distance from 
the rotation axis. 

The approximate alignment of the (warped) disk-like region of equilibria 
with the magnetic equatorial plane can be understood by rewriting Eq. Q in 
the form 



2 1- 



- 3(e,. • en) 



cos ip > ■ e,. = 



(16) 



If the star does not rotate near criticality, we can assume that up to some 
distance from the stellar surface the term in angular brackets dominates over 
the second term. In that limit, one solution of Eq. (|16|l represents the magnetic 
equator, em • e^ = 0. Consequently, for not too rapidly rotating stars, we can 
always expect an accumulation of matter in the magnetic equatorial plane near 
to the star. However, the deviation from a fiat disk increases for faster rotating 
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stars and is largest in the direction of the tilted dipole. Along these longitudes 
the disk is always bulged towards the rotational equator. 

3.4 Offset dipoles 

While we have treated so far only the case of a centered dipole, it appears as if 
this is more an exception than the rule for the type of stars we consider (e.g., 
Neiner 2002). A magnetic dipole with an offset vector a relative to the center 
of the star can be written as 

B„„.„, = -p- [3(ef • m)ef - m] , (17) 

where — (r — a)/f and f — |r — a|. The equilibrium positions and their 
stability are determined in full analogy to the case of centered dipoles. 

The qualitative behavior of the equilibrium positions for an offset dipole can 
be seen in Fig. O for the case -0 = 0.25 tt in comparison with the corresponding 
case with a centered dipole. The stellar parameters are the same as before. 
In the right panel of Fig. |2| the dipole has been moved by 0.5 i?* along the 
rotation axis (z-axis) . The shape of the equilibrium distribution is only slightly 
modified in comparison to the centered dipole. This is also the case for other 
realistic offsets (i.e., such that the dipole center remains within the star). Taking 
into account that we are mainly interested in the large scale distribution of 
matter ('^ 5 — 6 Tco) this result is not surprising since the field geometry on this 
scale is hardly infiuenced by dipole shifts ^ Itcq. Consequently, changes in 
the equilibrium distribution of circumstellar matter due to reasonable magnetic 
offsets are almost negligible. 

4 Discussion 

The origin of corotating circumstellar matter around hot stars with strong mag- 
netic fields is commonly explained with the models of Cassinelli et al. (2002) 
and Babel & Montmerle (1997) in terms of a radiatively driven stellar wind 
following the (approximately dipolar) magnetic field lines. Since the material 
streams from both hemispheres of the stellar surface and reaches supersonic 
speed, the model suggests that the plasma collides and forms a shock in the 
magnetic equator. The plasma cools through UV and X-ray emission and accu- 
mulates in an equatorial disk. While most of the theoretical calculations assume 
the special case of aligned magnetic and rotation axes, we have determined the 
equilibrium plasma distributions for the general case of an oblique dipolar ro- 
tator and also determined their stability properties. For aligned magnetic and 
rotation axes we confirm the accumulation of matter in a circumstellar disk. 
For increasing inclination angle the disk becomes warped and additional stable 
equilibrium positions appear. These are curved surfaces, roughly aligned with 
the axis of rotation, located above and below the magnetic equatorial plane. 
They grow and move towards the star for increasing values of ip. For the case 
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Figure 3: Equilibrium surfaces for inclination angle ^ = 0.25 tt. Left: Centered 
dipolc. right: dipole center displaced by along the rotation axis (^;-axis). 

Areas in light and dark gray indicate the stable and unstable equilibrium regions, 
respectively. 



of a decentered dipolc wc find that for reasonable offsets the stable equilibrium 
positions are changed only slightly compared with a centered dipole. We expect 
that, analogous to the formation of the equatorial disk, the newfound stable 
equilibrium positions lead to the local accumulation of plasma. Incoming ma- 
terial from the wind is cooled down by radiation behind shocks. Therefore, 
depending on the kinetic energy of the incoming plasma, we expect that the 
mass accumulation outside the disk appear in observations of continous and ro- 
tationally modulated IR, UV and X-ray emission, which is affected by the size 
of the emission region and its relative orientation to the line of sight. 

In fact, there are observations which indicate the existence of additional 
equilibrium positions outside the equatorial disk. For the rapidly rotating Be 
star u) Orionis, Neiner et al. (2003) found evidence for the presence of two coro- 
tating regions of material outside the plane of the disk by measuring variations 
in the peak intensities of emission lines. Such clouds were also reported by Short 
& Bolton (1994) for the magnetic He-strong star a Ori E. 

We are aware that our model is rather crude. Possible extensions include the 
dynamics of the radiatively driven wind and the corresponding cooling and accu- 
mulation processes. These should be incorporated into three-dimensional MHD 
simulations in order to understand the dynamics and distribution of plasma 
in rotating dipolar fields. First steps in this direction are the two-dimensional 
MHD simulations by Ud-Doula & Owocki (2002) (see also Ud-Doula & Owocki, 
2003). 
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5 Conclusion 



A simple consideration of the mechanical equilibrium of plasma in the magneto- 
sphere of a rotating magnetic star with relative inclination of the magnetic and 
rotation axes suggests the following distribution of matter accumulating, e.g., 
from a stellar wind: 

1. a disk- like structure roughly aligned with the magnetic equatorial plane 

2. two locations above and below the disk, coarsely aligned with the axis of 
rotation. 

The latter stable equilibrium regions are most prominent for large obliquity of 
the axes (perpendicular rotator). They provide hitherto unknown locations for 
the accumulation of circumstellar matter, which could explain some observa- 
tions of rapidly rotating magnetic stars, e.g. the close alignment of the new 
stable regions with the rotation axis could result in an observable rotational 
modulation of the disk radiation. 

We are grateful to Matthias D. Rempel and Francisco Frutos Alfaro for valuable 
discussions and to Coralie Neiner for pointing us to u) Orionis and aOri E. 
This research has made use of NASA 's Astrophysics Data System Abstract Ser- 
vice. 
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Figure 4: Equilibrium distribution for oblique rotators with various inclination 
angles (left panels: all equilibrium positions, right panels: only stable positions). 
Top: if) = 0.257r, middle: ^l) = 0. 37577, bottom: ij} = 0.457r. The rotation axis is 
always oriented along the 2;- axis. 



